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Abstract: The present investigation was carried out to study 
area, production and productivity trends and growth rates of 
maize (Zea mays) crop grown in Panchmahal region of Gujarat 
state, India for the period 1949-50 to 2007-08 based on 
parametric and nonparametric regression models. In parametric 
models different linear, non-linear and time-series models 
were employed. The statistically most suited parametric 
models were selected on the basis of adjusted R 2 , significant 
regression co-efficient and co-efficient of determination (R 2 ). 
Appropriate time-series models were fitted after judging the 
data for stationarity. The statistically appropriate model was 
selected on the basis of various goodness of fit criteria viz. 
Akaike’s Information Criterion, Bayesian Infonnation 
Criterion, RMSE, MAE , assumptions of normality and 
independence of residuals. In nonparametric regression 
optimum bandwidth was computed by cross-validation 
method. ‘Epanechnikov-kemeT was used as the weight 
function. Nonparametric estimates of underlying growth 
function were computed at each and every time point. Residual 
analysis was carried out to test the randomness. Relative 
growth rates of area, production and productivity were 
estimated based on the best fitted trend function. Linear model 
was found suitable to fit the trends in area and production of 
maize crop whereas for the productivity nonparametric 
regression without jump-point emerged as the best fitted trend 
function. The compound growth rate values obtained for the 
years (1949-50 to 2007-08) in area, production and 
productivity of the maize crop showed that the production had 
increased at a rate of 0.49 per cent per annum due to combined 
effect of increase in area and productivity at a rate of 0.30 and 
0.21 per cent per annum respectively. 

Key words: Adjusted R~, stationarity, akaike’s information 
criterion, bayesian information criterion, lijung and box test, 
cross validation, band width 

I. INTRODUCTION 

India is primarily an agriculture-based country and its 
economy largely depends on agriculture. About 25% of our 
country’s Gross Domestic Product (GDP) comes from 
agricultural sector. Nearly 74 % of the country’s population 
lives in villages and depends on agriculture. Therefore, 
country’s development largely depends upon the development 
of agriculture (Nath, 2008). The statistical information on crop 
area, production and productivity form the backbone of 
agricultural statistical system. Regional data analysis is 
extremely vital since, it forms the basis for economic and 
policy planning by the state and central governments. It is easy 
to formulate and initiate appropriate policy measures if the 
data with regard to the trend of production is obtained and 
analysed in advance. Growth rate analyses are widely 
employed to study the long-term trends in various agricultural 
crops (Panse, 1964). 


The growth rates of different crops are estimated 
mostly through the parametric models by assuming the linear 
or exponential functional forms. A number research workers 
(Panse, 1964; Dey, 1975; Reddy, 1978; Narain et al., 1982; 
Patel et al., 1986, Kumar and Rosegrant, 1994; Kumar, 1997; 
Borthakur and Bhattacharya, 1998; Joshi and Saxena, 2002; 
Sharma, 2005; Patil et al., 2009) have used parametric models, 
to estimate growth rates, which are currently being used by the 
planners or policy makers of the country. However, the data 
may not be following these linear or exponential models or 
may require fitting of higher degree polynomials or non-linear 
models. Further this models lack the econometric 
consideration i.e., normality and randomness of residuals. 
Under these circumstances it becomes imperative to take 
recourse to nonparametric regression approach, which is based 
on fewer assumptions. 

The objective of the present study is to develop an 
appropriate statistical model to fit the trends and to calculate 
growth rates in area, production and productivity of maize crop 
grown in Panchmahal region of Gujarat state based on both 
parametric (Linear, non-linear and time-series) and 
nonparametric regression models. 

II. MATERIALS AND METHODS 

To achieve the stipulated objectives, the present study 
had been carried out on the basis of time-series data pertaining 
to the period 1949-50 to 2007-08 had been collected from 
various publication (Margdarshika, published yearly by 
Directorate of Agriculture) of Gujarat government 
(Anonymous, 2009). 

In parametric models different linear (Montgomery et 
al., 2003), non-linear (Ratkowsky, 1990; Bard, 1974; Draper 
and Smith, 1998) and Auto-Regressive Integrated Moving 
Average (ARIMA) time-series models (Box et al.. 1976) were 
employed. The statistically most suited parametric models 
were selected on the basis of adjusted R 2 , significant regression 
co-efficient, co-efficient of determination (R 2 ), Root Mean 
Square Error (RMSE), Mean Absolute Error (MAE) and 
assumptions of residuals (normality and randomness). 

Appropriate ARIMA models were fitted after judging 
the time-series data for stationarity based on visual inspection, 
auto-correlation function and partial auto-correlation function. 
The auto-correlation upto fifteen lags were worked out. The 
statistically most appropriate time-series model was selected 
based on various goodness of fit criteria viz. Akaike’s 
Information Criterion (AIC), Bayesian Information Criterion 
(BIC), RMSE, MAE and assumptions of residuals (Shaprio- 
Wilk test for normality and Lijung & Box test for 
randomness). 

In nonparametric regression (Hardle, 1990), the first 
step involved estimation of optimum bandwidth and was 
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computed by cross-validation method. ‘Epanechnikov-kernel’ 
was used as the weight function. Nonparametric estimates of 
underlying growth function were computed at each time point. 
Residual analysis was carried out to test the randomness. A 
relative growth rate was calculated based on best fitted model. 

III. RESULTS AND DISCUSSION 

Different linear, non-linear and nonparametric 
regression models were employed to study the trends in area, 
production and productivity of maize crop. Data obtained for 
Panchmahal district of middle Gujarat were employed. The 
characteristics of fitted linear, non-linear (Tables 1, 2. and 3) 
and time-series (Table 4) models are presented. The findings 
are discussed in sequence as under. 

Trends in area, based on linear and non-linear models: 

The data presented in Table 1 for area under the 
maize crop revealed that among the linear and non-linear 
models fitted for area under the maize crop, the maximum 
adjusted R 2 of 91 per cent was observed in case of third degree 
polynomial model and since the run test value was found to be 
significant, this model failed to fulfill the model selection 
criteria. The quadratic, exponential and Gompertz models also 
failed to fulfill the model selection criteria. The linear model 
with the maximum of 85 per cent adjusted R 2 with the 
comparatively lower values of RMSE (120.60) and MAE 
(87.86) was found suitable to fit the trends in area under the 
cotton crop. 

Y= 1069.76** + 17.19** X 

(R 2 = 85** %) 

Trends in area, based on time-series models: 

In ARIMA time-series methodology the auto¬ 
correlation upto fourteen lags was worked out. Since the 
computed auto-correlations y k values did not tail off towards 
zero, the original series was found to be non-stationary. The 
non-stationarity was also confirmed by examining the 
realization visually. It was found that the mean and variance 
were changing over the time. However, the stationarity was 
achieved by differencing two times i.e., d=2. The pattern of 
auto-correlations y k showed damped sine-wave and significant 
partial auto-correlations (|)kk at second and third lags. This 
suggested consideration of ARIMA(2,2,0) and ARIMA(3,2,0) 
as the candidate models and the results are given in Table 4. 
Since the Shapiro-Wilks test was significant in both the 
models, these models have failed to fulfill the model selection 
criteria. Hence none of the ARIMA families’ of time-series 
models were found suitable to fit the trend in area under the 
maize crop. 

Trends in area, based on nonparametric regression model: 

Nonparametric regression could not be fitted to the 
area under the maize crop due to over smoothing. As per the 
model selection criteria linear model was selected as the best 
fitted model to fit the trends in area under the maize crop. The 
graph of the fitted trends in maize crop using linear model is 
depicted in the Fig. 1. 

Trends in production, based on linear and non-linear 
models: 

The data presented in Table 2 for production of maize 
crop revealed that among the linear and non-linear models 
fitted, the maximum adjusted R 2 of 26 per cent was observed 
in case of linear model and this model showed comparatively 
lower values of root mean square (767.52) and mean absolute 


(613.20) errors in comparison to that of other linear and non¬ 
linear models. Since all the partial regression co-efficient was 
highly significant and Shapiro-Wilks test and run test values 
were non-significant, this model fulfilled model selection 
criteria. The following model was found suitable to fit the 
trend in production of maize crop. 

Y = 869.85“ + 27.54** X 
(R 2 = 27** %) 

Trends in production, based on time-series models: 

In ARIMA time-series methodology the auto¬ 
correlation upto fifteen lags was worked out. Since the 
computed auto-correlations yk values did not tail off towards 
zero, the original series was found to be non-stationary. The 
non-stationarity was also confirmed by examining the 
realization visually. It was found that the mean and variance 
were changing over the time. However, the stationarity was 
achieved by differencing one time i.e., d=l. The pattern of 
auto-correlations y k showed damped sine-wave and significant 
partial auto-correlations 4>kk at third, fourth and fifth lags. This 
suggested consideration of ARIMA(3,1,0), ARIMA(4,1,0) and 
ARIMA(5,1,0) as the candidate models and the results are 
given in Table 4. But the values of AIC (954.14), BIC 
(966.51), RMSE (851.72) and MAE (634.31) were minimum 
in ARIMA(5,1,0) model in comparison to that of other 
candidate model. Also the model ARIMA(5,1,0) fulfilled the 
model selection criteria and hence the model ARIMA(5,1,0) 
was found suitable to fit the trend in production of maize crop 
among the ARIMA families’ of time-series models. 

Trends in production, based on nonparametric regression 
model: 

Nonparametric regression model could not be fitted in 
case of production of maize crop due to over smoothing. Since 
the RMSE (767.52) and MAE (613.20) values in linear model 
were minimum in comparison to those in ARIMA (5,1,0), the 
linear model was found suitable to fit the trend in production 
of maize crop. The graph of the fitted trend for production of 
maize crop using linear model is depicted in the Fig.2. 

Trends in productivity, based on linear and non-linear 
models: 

The data presented in Table 3 for productivity of 
maize crop revealed that among the linear and non-linear 
models fitted, the adjusted R 2 ranges from 01 to 05 per cent. 
These values were too low to fit the trends in productivity of 
maize crop. This might be due to very larger variability in 
yield of maize crop. None of the linear and non-linear model 
was found suitable to fit the trends in productivity of maize 
crop. 

Trends in productivity, based on time-series models: 

None of the ARIMA families’ time-series models was 
found suitable to fit the productivity trend in maize crop. 

Trends in productivity, based on nonparametric regression 
model: 

For the productivity data of maize crop the optimum 
bandwidth was computed as 0.07 using the cross-validation 
method. Nonparametric estimates of underlying growth 
function were computed at each and every point. Residual 
analysis showed that the assumptions of independence of 
errors were not violated at 5% level of significance. The root 
mean square and mean absolute errors values were found to be 
425.80 and 360.35 respectively. The nonparametric regression 
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model was selected as the best fitted trend function for the 
productivity trend of the maize crop. No significant jump-point 
was observed in the productivity of maize crop. The graph of 
the fitted trend for the productivity of maize crop using the 
nonparametric regression is depicted in the Fig.3. 

Discussion in area, production and productivity of paddy 
crop: 

Yadav and Das (1990) and Bera et al., (2002) used 
exponential model to study the trends in area, production and 
productivity of maize in Assam for the periods from 1966 - 67 
to 1985 - 86 and in major maize growing districts of West 
Bengal respectively. In the present study exponential model 
could not be emerged as the best fitted model because of lack 
of assumptions of residuals. The linear model emerged as the 
best fitted trend model for area and production of maize crop. 

Nonparametric regression was used to fit the trends in 
productivity of maize crop 

Growth rates in area, production and productivity of paddy 
crop: 

The compound growth rate values obtained for the 
years (1949-50 to 2007-08) in area, production and 
productivity of the maize crop showed that the production had 
increased at a rate of 1.82 per cent per annum which might be 
due to combined effect of increase in area and productivity at a 
rate of 1.12 and 0.77 per cent per annum respectively. Growth 
rate of area, production and productivity of maize crop were 
depicted in Fig.4. 



Figure 4: Compound Growth Rates of area, production and 
productivity of maize crop 
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Table 1: Characteristics of fitted linear and non-linear models for area of maize crop (Panchmahal) 


Model 

Regression Coefficients 

Goodness of fit 

A 

B 

C 

D 

R 2 / 

Adj.R 2 

Shapiro 

-Wilks 

test 

Run 

Test 

RMSE 

MAE 

Linear 

1069.76” 

(32.36) 

17.19" 

(0.938) 

- 

" 

0.85” 

[0.85] 

0.087 

0.059 

120.60 

87.86 

Quadratic 

1213.59" 

(43.13) 

3.04 

(3.32) 

0.2358" 

(0.0536) 


0.89" 

[0.89] 

0.006 

0.088 

103.96 

74.89 

Cubic 

1085.15" 

(54.25) 

27.70" 

(7.76) 

-0.7831” 

(0.2993) 

0.0113" 

(0.0033) 

0.91" 

[0.91] 

0.237 

0.006 

94.26 

71.24 

Exponential 

1128.89" 

(20.32) 

0.0107" 

(0.000522) 

- 


0.88" 

[0.88] 

0.016 

0.049 

110.47 

79.87 

Monomolecular 

52572.40’ 

(1099933.99) 

51503.84' 

(1099895.94) 

0.000337' 

(0.00726) 

- 

0.85' 

[0.84] 

0.086 

0.049 

120.98 

88.08 

Gompertz 

Relation 

672712.92' 

(14618521.75) 

1.86’ 

(3.39) 

0.0018' 

(0.0066) 

- 

0.88" 

[0.88] 

0.041 

0.013 

111.26 

81.95 


Table 2: Characteristics of fitted linear and non-linear models for production of maize crop (Panchmahal) 


Model 

Regression Coefficients 

Goodness of fit 

A 

B 

C 

D 

R 2 / 

Adj.R 2 

Shapiro 
- Wilks 

test 

Run 

Test 

RMSE 

MAE 

Linear 

869.85" 

(205.93) 

27.54" 

(5.97) 

- 

- 

0.27" 

[0.26] 

0.879 

0.894 

767.52 

613.20 

Quadratic 

931.34" 

(318.24) 

21.49 

(24.47) 

0.1008 

(0.3954) 

- 

0.27" 

[0.25] 

0.846 

0.897 

767.07 

612.29 

Cubic 

767.15 

(440.30) 

53.01 

(63.02) 

-1.2016 

(2.4296) 

0.0145 

(0.0266) 

0.28" 

[0.24] 

0.773 

0.897 

765.02 

610.67 

Exponential 

891.49" 

(136.01) 

0.0158" 

(0.0044) 

- 

- 

0.18" 

[0.17] 

0.987 

0.897 

799.18 

652.41 

Monomolecular 

39924.91' 

(1529527.82) 

39060.56* 

(1529286.32) 

0.00071 

(0.0287) 

- 

0.27* 

[0.25] 

0.883 

0.894 

768.00 

614.00 

Gompertz 

Relation 

21941.16' 

(257639.72) 

1.1419* 

(3.67) 

0.0064* 

(0.0298) 

- 

0.27* 

[0.24] 

0.841 

0.359 

767.12 

611.92 

Hoerl 

794.67* 

(414.62) 

1.012* 

(0.0010) 

0.1171* 

(0.2410) 

- 

0.28* 

[0.25] 

0.840 

0.897 

765.56 

611.31 

Logistic 

6852.00' 

(27419.49) 

6.085’ 

(26.99) 

0.0220* 

(0.0311) 

- 

0.27* 

[0.24] 

0.840 

0.359 

767.18 

611.81 


* Significant at 5% level ** Significant at 1% level 

Values in brackets ( ) indicate standard errors Values in square brackets [ ] indicate 

Adjusted R 2 


Table 3: Characteristics of fitted linear and non-linear models for productivity of maize crop (Panchmahal) 


Model 

Regression Coefficients 

Goodness of fit 

A 

B 

C 

D 

R 2 / 

Adj.R 2 

Shapiro 
- Wilks 

test 

Run Test 

RMSE 

MAE 

Linear 

856.18" 

(115.67) 

6.25 

(3.35) 

- 

" 

0.06 

[0.04] 

0.115 

0.695 

431.12 

361.77 

Quadratic 

744.03" 

(177.76) 

17.27 

(13.67) 

-0.1839 

(0.2208) 

- 

0.07 

[0.04] 

0.111 

0.897 

428.47 

361.49 

Cubic 

658.22" 

(246.03) 

33.75 

(35.21) 

-0.8646 

(1.3576) 

0.0076 

(0.0149) 

0.07 

[0.02] 

0.109 

0.897 

427.47 

359.70 

Exponential 

789.71" 

(116.06) 

0.0051 

(0.0043) 

- 

- 

0.02 

[0.01] 

0.150 

0.695 

448.39 

386.93 

Monomolecular 

1110.00* 

(75.56) 

667.85* 

(408.03) 

0.156* 

(0.146) 

- 

0.08* 

[0.05] 

0.122 

0.894 

425.51 

358.50 

Gompertz 

Relation 

1105.15* 

(70.89) 

-0.1087* 

(0.8614) 

0.2039* 

(0.1897) 


0.08* 

[0.051 

0.143 

0.894 

422.39 

358.28 

Logistic 

1102.25* 

1.37* 

0.2516* 

- 

0.08* 

0.161 

0.510 

425.30 

358.23 
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(68.51) 

(1.58) 

(0.2367) 


10.051 





Fioerl 

594.71* 

(231.59) 

0.9965* 

(0.0086) 

0.2106’ 

(0.1903) 

- 

0.08’ 

[0.05] 

0.098 

0.897 

426.02 

359.23 

Morgan- 

Mercer-Flodin 

541.80* 

(415.76) 

58.50* 

(477.44) 

1110.48* 

(92.43) 

2.37* 

(3.83) 

0.09* 

[0.04] 

0.157 

0.894 

424.63 

357.32 


Table 4: Characteristics of fitted time-series models for area, production and productivity of Maize (Panchmahal) 


Aspects 

ARIMA 

(P,d,q) 

Constant 

Auto-regressive Co-efficient 

Goodness of Fit 

<t>i 

<l>2 

4*3 

4*4 

<l>5 

4*6 

AIC / 

BIC 

Shapiro 

-Wilks 

test 

Box - 
Ljung 

RMSE/ 

MAE 

Area 

(2,2,0) 

-0.868 

-1.051 

-0.414 

- 

- 

- 

- 

716.80 

0.001 

8.48 

140.39 



(6.840) 







/ 


[21.02] 

/ 97.41 




(0.129) 

(0.124) 





722.93 





(3,2,0) 

-0.918 

-1.108 

-0.592 

-0.189 

- 

- 

- 

717.14 

0.001 

9.32 

140.07 



(5.821) 



(0.137) 




/ 


[19.67] 

/ 94.83 




(0.139) 

(0.188) 





725.32 




Production 

(3,1,0) 

21.333 

-0.669 

-0.376 

-0.352 

- 

- 

- 

959.83 

0.411 

10.60 

914.09 



(50.832) 







/ 


[19.67] 

687.80 




(0.128) 

(0.150) 

(0.136) 




968.07 





(4,1,0) 

22.925 

-0.732 

-0.471 

-0.469 

-0.212 

- 

- 

959.65 

0.137 

9.94 

904.55 



(42.011) 

** 



(0.149) 



/ 


[18.30] 

/ 




(0.134) 

(0.162) 

(0.158) 




969.96 



685.52 


(5,1,0) 

28.103 

-0.804 

-0.683 

-0.689 

-0.499 

-0.399 

- 

954.14 

0.501 

5.65 

851.72 



(28.263) 

** 

** 

** 

** 

** 


/ 


[16.91] 

/ 




(0.127) 

(0.170) 

(0.168) 

(0.176) 

(0.152) 


966.51 



634.31 

Productivity 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 

- 


** estimated t - values are greater than or equal to 2 Values in the brackets () are corresponding 

standard errors 

Values in the square brackets [ ] indicate critical values for Chi - square statistic at 5 % level of significance 



1 3 5 7 9111315171Q212725X724313335373941434547495153555751 

X ™.. 1 


Figure 1: Trends in area of maize crop based on first degree polynomial model 



Figure 2: Trends in production of maize crop based on first degree polynomial 
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Figure 3: Trends in Productivity of Maize Crop Based On Nonparametric Regression 


CONCLUSION 

Linear model was found suitable to fit the trends in 
area and production of maize crop whereas for the productivity 
nonparametric regression without jump-point emerged as the 
best fitted trend function. The compound growth rate values 
obtained for the years (1949-50 to 2007-08) in area, production 
and productivity of the maize crop showed that the production 
had increased at a rate of 0.49 per cent per annum due to 
combined effect of increase in area and productivity at a rate of 
0.30 and 0.21 per cent per annum respectively. 
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